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摘 要 界面 声 子 热 输 运 现象 广泛 地 存在 于 纳米 热电 装置 和 微 纳 电 子 器 件 中 , 但 其 物理 机 制 和 介 观 模型 仍 不 完善 . 本 文 
将 考虑 真实 色散 关系 的 漫 射 失 配 (DMM) 界面 模型 引入 到 高 效 、 低 噪 的 动力 型 声 子 Monte Carlo(MC) 方法 , 建立 了 一 个 
频谱 穿 透 系数 的 界面 声 子 输 运 数值 模型 . 采用 新 的 数值 模型 研究 了 室温 下 Si/Al 薄膜 中 的 界面 声 子 输 运 , 可 准确 捕捉 法 向 


输 运 和 面向 输 运 的 非 平衡 效应 。 


关键 词 


中 图 分 类 号 : TK124 文献 标识 码 : A 


声 子 输 运 ， Monte Carlo 模拟 ; 界面 热 输 运 ; 微 纳 尺度 传 热 
文章 编号 : 0253-231X(2017)04-0841-06 


Monte Carlo Simulation of Interfacial Phonon Transport Through 
Bilayer Thin Fil 


RAN Xin 


GUO Yang-Yu 


(Key Laboratory for Thermal Science and Power Engineering of 


and CNMM, Tsinghua U. 


sit, 


AAA 
G Mo-Ran 


istry of Education, Department of Engineering Mechanics 
, Beijing 100084, China) 


Abstract Interfacial phonon transport itt exists in nano thermoelectric and nano electronic 
devices, but the mechanism and a aa remain still incomplete. The present work estab- 


lishes a numerical model for interfi 
(DMM) with exact phonon dis 


to a kinetic-type Monte Carlo 


onon transport by introducing the diffuse mismatch model 


(MC) scheme. Then the interfacial 


phonon transport through Si (A Dilayer thin film at room temperature is studied. The new numeri- 


ell the non-equilibrium effects in both cross-plane and in-plane transport. 


cal model is able to ae 
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0 引 


近 几 十 年 来 ， 随 着 微 纳 加 工 制造 技术 和 微 纳米 
科技 的 发 展 , 微 纳 尺度 热 输 运 成 为 人 们 的 研究 热点 
在 微 纳 系统 中 的 热 输 运 ， 界 面 效应 占 主导 作 寺 兴 适 
用 于 体 材料 的 Fourier 导热 定律 不 再 成 立 和 出， 因而 
耻 需 发 展 新 的 理论 和 模型 . 在 传统 的 半导体 器 件 中 ， 
声 子 和 电子 是 主要 的 热 载 子 (Heat carrier)， 而 本 文 
主要 关注 声 子 热 输 运 。 微 纳 尺 度 热 输 运 领域 主要 由 
田 长 霖 等 人 在 上 世纪 八 九 十 年 代 开创 P, ZIJL 
年 的 发 展 ， 已 有 十 足 的 进展 ， 但 仍 存 在 很 多 尚未 解 
决 的 难题 ,例如 微 纳 热电 材料 的 能 量 转换 效率 优化 
和 微 纳 电子 器 件 散热 . 这 两 种 难题 中 都 存在 共同 的 
界面 声 子 热 输 运 过 程 . 

前 人 在 界面 声 子 输 运 方面 进行 了 大 量 的 理论 研 
究 ， 提 出 的 经 典 界 面 输 运 模型 主要 包括 : 声学 失 配 
模型 (Acoustic mismatch model, AMM) 中 和 漫 射 失 
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oo micro- and nano- 


Æ] (Diffuse mismatch model, DMM) 4], AMM 
模型 通常 只 适用 于 极 低温 情形 ， 认 为 声 子 以 连续 波 
的 形式 在 介 电 固体 中 传播 ， 并 在 光滑 界面 发 生 镜 面 
反射 和 镜面 透射 ， 穿 透 系数 (透射 率 ) 可 由 界面 两 侧 
材料 的 声学 阻抗 计算 获得 ， 反 射 角 和 透射 角 由 光学 
中 的 Snell 定律 计算 获得 加 ;而 DMM 模型 适用 于 
更 高 温度 的 情形 ， 认 为 声 子 在 界面 发 生 漫 反 射 和 漫 
透射 ， 穿 透 系数 主要 取决 于 界面 两 侧 材 料 的 声 子 色 
BOCA |), 前 人 的 实验 研究 间接 表明 了 声 子 在 界面 
处 存在 镜 射 和 漫 射 行为 ,因此 AMM 模型 和 DMM 
模型 均 无 法 单独 地 准确 描述 所 有 的 界面 声 子 输 运 过 
程 回 .Chen (9 在 1998 年 提出 一 个 基于 Debye 近似 
的 灰色 混合 模型 ,假设 声 子 在 界面 发 生 部 分 镜 射 和 
部 分 漫 射 , 但 镜 参 数 ( 镜 射 的 比例 ) 是 一 个 经 验 值 , 通 
稼 需要 通过 与 实验 结果 的 对 比 获 得 ， 声 子 界面 输 运 
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行为 的 本 质 仍 有 等 进一步 的 揭示 。 由 于 人 们 感 兴 
的 微 纳 尺度 热 输 运 过 程 通 常 发 生 在 室温 及 更 高 温度 
F, DMM 模型 为 目前 较为 广泛 认可 和 使 用 的 模型 
后 人 在 DMM 模型 的 基础 上 也 提出 了 很 多 的 修正 或 
者 改进 『- 世 0， 虽 然 无 法 准确 地 描述 所 有 的 声 子 界面 
输 运 , 但 更 近 一 步 地 解释 了 实验 结果 21, 

结合 上 述 界面 输 运 模型 ， 前 人 通过 求解 声 子 玻 
尔 兹 曼 方 程 对 微 纳 尺度 界面 声 子 输 运 进行 了 很 多 的 
研究 . 求解 声 子 玻 尔 效 曼 方程 的 方法 主要 包括 确定 
型 (有 限 差 分 法 FDM"4) 等 ) 和 随机 型 (Monte Carlo, 
MC [5-171) 两 种 . MC 方法 避免 了 直接 求解 高 维 的 声 
子 玻 尔 兹 曼 方 程 ， 且 对 界面 边界 的 处 理 简 单 并 具有 
清晰 的 物理 图 像 ， 因 而 是 复杂 界面 (多 层 薄 膜 结构 、 
纳米 复合 材料 和 纳米 唱 等 ) 中 声 子 输 运 研究 的 良好 
选择 .Jeng $US 最 早 使 用 声 子 MC 方法 计算 了 纳 
米 颗 粒 复合 材料 的 等 效 热 导 率 ,采用 了 基于 Debye 
近似 的 灰色 DMM 模型 ;后 来 Huang 等 65] 在 声 子 
MC 中 考虑 了 Chen 提出 的 部 分 镜 射 和 部 分 漫 射 的 灰 


使 用 新 发 展 的 能 量 偏差 声 子 MC 方法 模拟 了 

反射 实验 中 Si/Al 界面 的 多 尺度 声 子 输 运 ; 采 
用 了 Minnich 等 07) 提出 的 修正 DM 。 总 之 ， 
前 人 的 介 观 模拟 工作 在 界面 边界 处 理 时 均 只 考虑 了 


色 混 合 界面 模型 ; 最 近 Peraud 和 ee 


常数 型 声 子 界面 穿 透 系数 ,或 从 一 个 方向 穿 
过 界面 时 的 穿 透 系 数 是 一 个 党 数 . 但 最 近 的 实验 测 


量 8) 和 微观 理论 计算 结果 [9 均 表 明 , 声 子 的 界面 
穿 透 系数 具有 高 度 的 频谱 特性 ， 通常 而 言 ， 高 频 声 
子 在 界面 处 容易 被 反射 ， 而 低频 声 子 则 较 容 易 通 过 
界面 . 前 人 的 模型 研究 也 表明 考虑 真实 色散 关系 和 
频谱 穿 透 系数 的 DMM 模型 与 Debye 灰色 DMM 模 


型 相 比 更 符合 实际 ri， 计算 获得 的 界面 热 导 与 实 。 


验 结果 更 接近 Da 然而 ,目前 仍然 缺少 一 个 考虑 
谱 穿 透 系数 的 界面 声 子 输 运 介 观 MC BRON a 
此 , 本 文 将 基于 考虑 真实 色散 关系 的 DMM 模型 , 建 
立 一 个 考虑 频谱 穿 透 系数 的 双 层 薄膜 界面 声 子 输 运 
MC 数值 模型 ， 


1 物理 模型 和 数值 方法 


1.1 物理 模型 

本 文 以 室温 下 硅 / 铝 薄膜 界面 声 子 输 运 为 研究 对 
RR, 考虑 了 法 向 输 运 和 面向 输 运 两 种 情形 , 分 别 如 图 
1(a)、 图 1(b) Bram. 对 于 法 向 输 运 , 系统 处 于 均匀 的 
初始 温度 299 K， 突 然 系统 的 左边 界 与 301 K 的 恒 
温 热 源 接触 , 右边 界 与 299 K 的 恒温 热源 接触 .y TF 
向 上 尺度 无 限 大 ， 以 保证 一 维 的 法 向 输 运 。 对 于 面 
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向 输 运 , 在 x 方向 上 施加 一 个 恒定 的 均匀 温度 梯度 ， 
上 下 边界 绝热 


无 限 大 
(a) 法 向 输 运 
(a) Cross-plane transport 


绝热 


(b) 面向 输 运 
(b) In-plane transport 
图 1 硅 / 铝 薄膜 界面 声 子 输 运 物理 模型 示意 图 
Fig. 1 Physical model of Si/Al thin flm interfacial phonon 


transport 


本 文采 用 硅 、 铝 体 材料 的 真实 的 声 子 色散 关 
系 Po, 忽略 光学 支 对 热 输 运 的 贡献 ， 且 使 用 四 次 
多 项 式 拟 合 表达 式 ， 即 : 


WwW (k) = Agk* + Ask? + Aok? + Aik (1) 


AH, w 为 声 子 频率 , k HEATH, Ai, Ao. As 
和 Ay 为 拟 合 参数 。 对 于 硅 体 材 料 ， 拟 合 参数 分 别 
为 : 横 声 学 声 子 (TA): A1=5.511x103 msl, A= 
一 1.169x10-8 m?-s-!, Ag=—4.957x107)" m3.s-1, 

44=2.432x10-27 m4.s-1， 纵 声学 声 子 (LA): A= 
8.192x103 m-s7!, Ag=—1.140x1077 m2.s-1，43 = 
一 2.612 x 10718 m3 . s-1， As=—5.645x10-29 m4. 


4 期 


s 1 [29。 对 于 铝 体 材料 ， 考 虑 真实 的 色散 关系 PN, 
我 们 同样 采用 四 次 多 项 式 拟 合 ,我 们 拟 合 得 到 的 无 
基 纲 的 色散 关系 表达 式 为 : 

w(K) = AgK*+A3K3+ AK? +AiK (2) 


式 中 , K = k/kmax, kmax=27/a, a 为 铝 的 唱 格 常 数 ， 
拟 合 参数 分 别 为 : 横 声 学 声 子 (TA): 41=5.147x1073 
rad/s, 42=9.578x1012 rad/s, A3=—3.204x10!% 
rad/s, A4=6.750x10! rad/s; MEZAT (LA): 
A1=9.825x1013 rad/s, A2=2.226x10!3 rad/s, A3= 
—1.216x1014 rad/s, Ay=6.208x 101% rad/s. 

FED BLA ee eT SEF Matthiessen 定 
W, Aer tatty, 其 中 非 简 谐 三 声 子 过 程 的 
HRT ARIAT Pl. 


7} = ATW’ exp(—B/T) (3) 


AH, A=1.4x107!9 s/K, B=152 K。 缺 陷 散 射 过 程 
EY a BRAY Ta EIA PO; 


“t= Cut 


AF, C=1.32x10-* s3, FARR AE 
AMM: 7 = 1x10-11 s, 声 子 的 界面 输 i 考虑 真 
: 
> 
or (w) = er _ (5) 


2 lkp YP + > Thp.2 ©) 


AH, 下 标 1 和 2 代表 材料 , p RRA FX. 根据 漫 
射 假设 , A: az 1 (w) = 1- ar (w). 
1.2 声 子 蒙特 卡 洛 方法 

本 文采 用 Peraud 和 
发 展 的 动力 型 Monte Carlo 方法 [22， J 
用 a E.G 


方程 : 
o (e zA) 十 vlw, p) a V(e en )= = 
(ee = ft.) 一 (c -er J (6) 


T (w, p, T) 


AF, e 为 声 子 的 能 量 分 布 ，vws 为 声 子 群 速度 , w 为 
ETAMZ, p ARP, Toy 为 假设 的 参考 平衡 温 
BE, en, 为 在 Teg 下 的 平衡 能 量 分 布 ， eloc 为 在 当地 
伪 温 度 Toc 的 平衡 能 量 分 布 , 7 WETERE. 基 
于 能 基 偏 差 的 MC 方法 只 模拟 导热 过 程 偏离 参考 平 
衡 分 布 的 声 子 分 布 , 因而 统计 噪音 小 ; 在 进一步 的 小 
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温差 假定 下 , 动力 型 MC 将 偏差 能 量 的 分 布 做 线性 
化 近似 , 归 一 化 后 为 与 当地 温度 无 关 的 分 布 , 因而 无 
需 在 中 间 步 又 中 对 宏观 变量 (温度 和 伪 温 度 ) 进行 统 
计 来 获得 过 程 中 声 子 的 平衡 分 布 和 非 平 衡 分 布 ， 只 
要 知道 所 需 发 射 的 声 子 的 总 数目 ， 就 可 以 一 个 一 个 
PRES, 相 比 通常 的 MC E 进一步 简化 了 数值 算法 ， 
提高 了 计算 效率 。 
动力 型 声 子 MC 方法 的 主要 流程 依次 包括 初始 
化 、 迁 移 、 散 射 、 边 界 处 理 和 宏观 量 统计 . 首先 是 初 
台 化 步 ， 给 定 模拟 粒子 如 下 的 平 衔 分 布 : 
der. 
ar (7) 
To 为 系统 初始 温度 , Dw, p) 为 声 子 态 密度 , 粒子 以 
平衡 分 布 (7) 从 边界 或 系统 内 部 一 个 一 个 地 依次 发 
aT ORR mae 模拟 粒子 的 迁移 位 移 是 VAt, 
—T-In(R), t 是 模拟 时 间 , 7 是 材料 的 声 
Ro R Ge 0 到 1 的 随机 数 ， 迁移 
藻 粒 子 遇 到 边界 , 则 进行 边界 处 理 , 否则 对 粒子 进 
站 行 散射 处 理 ， 散射 后 的 粒子 服从 的 非 平衡 分 布 为 : 
D(w,p) 
T(w, P, Teq) dT 


一 次 散射 结束 后 , 粒子 继续 进行 漂移 和 下 一 次 散射 ， 
如 此 循环 往复 ， 直 至 时 间 到 达 设 定 的 计算 终止 时 间 
或 者 粒子 离开 模拟 系统 边界 。 最 后 对 时 间 和 空间 进 
行 离散 化 ， 并 在 合作 个 时 间 节 点 上 对 声 子 在 空间 网 
piles 总 能 量 流 进行 统计 以 获得 温度 和 热 
ane en 最 后 对 边界 处 理 进行 说 明 ， 本 文 涉 及 
` 理 包括 绝热 边界 、 等 温 边界 、 恒 定 温度 梯 

和 界面 边界 . 绝热 边界 采用 漫 反 射 边 界 处 理 ， 


(To = Tea) D(w, p) 


der, 


(Toc ~ Teq) 


(8) 


in mon anaes 


与 入 射 粒子 的 速度 方向 无 关 , 但 是 粒子 的 频率 和 波 
数 保持 不 变 。 等 温 边 界 要 求 与 边界 磁 撞 的 粒子 全 吸 
收 , 发 射 粒子 是 漫 发 射 , 即 方向 随机 给 定 , 满足 平衡 
分 布 (7). 面向 输 运 中 的 恒定 温度 梯度 边界 采用 热流 
周期 性 边界 处 理 ， 从 一 个 边界 离开 系统 的 粒子 从 另 
外 一 个 边界 重新 进入 系统 ， 除 位 置 之 外 的 其 它 粒子 
信息 都 不 改变 . 界面 边界 采用 漫 射 处 理 ， 随 机 生成 
一 个 0 到 1 之 间 的 数 ， 并 与 穿 透 系数 比较 ， 若 小 于 
穿 透 系数 , 则 漫 穿 透 ; 若 大 于 穿 透 系 数 ， 则 漫 反 射 。 
2 程序 验证 

采用 前 面 介绍 的 动力 型 MC 方法 , 模拟 了 室温 
(T=300 K) 下 单 层 硅 薄膜 中 声 子 输 运 过 程 ， 考 虑 了 
法 向 输 运 和 面向 输 运 两 种 情形 , 分 别 与 解析 结果 进 
行 对 比 , 验证 程序 的 正确 性 . 
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对 于 硅 薄 膜 中 的 法 向 声 子 输 运 ,为 了 方便 与 已 
有 法 向 输 运 的 结果 PS) 进行 对 比 , 采用 声 子 灰色 模型 
并 作 Debye 近似 , 单 唱 硅 体 材料 的 热 物性 参数 如 表 
1 所 示 多 ,法 向 Kn 数 ( 声 子平 均 自由 程 与 薄膜 厚度 
之 比 ) 依次 为 0.035、0.35 和 3.5, 统计 宏观 量 的 时 间 
间隔 分 别 为 5x10-13 s、5x10-13s 和 3x10-14s, 网 
格 数目 是 400、40 和 40, 计算 总 时 间 分 别 是 6x10-8 
s, 1.2x107-° s 和 3x10-1 s, 偏差 能 量 粒 子 数目 是 
6000 万 、3000 万 和 4000 万 . 法 向 声 子 输 运 的 模拟 结 
果 如 图 2 所 示 . 图 2(a) RH, 稳 态 温度 分 布 均 近 似 
为 线性 , 但 随 着 Kn 的 增 大 , 在 边界 处 的 温度 跳跃 越 
明显 , 法 向 非 平衡 效应 越 强 , 热 阻 越 大 , 因而 法 向 有 
效 热 导 率 越 小 ， 如 图 2(b) Pras. 当前 的 温度 分 布 和 
法 向 有 效 热 导 率 的 MC 模拟 结果 均 与 文献 中 已 有 的 
数值 和 理论 结果 吻合 很 好 . 


表 1 硅 晶 体 在 300 K 的 热 物理 性 质 


Table 1 Thermophysical properties of Si TY 


at 300 K 
Cv/(J/(m3K)) T, 
1.18x 10° 


1 


l/nm 
41.79 


Ug /m-s7 


o MC simulation 
(b) 
10! 10° 107! 107? 107° 
Kn 
图 2 室温 下 硅 薄 膜 中 的 稳 态 法 向 声 子 输 运 的 MC 模拟 结 
Fe: (a) 不 同 Kn 数 下 的 无 量 纲 温度 分 布 , 符号 是 当前 的 MC 数 
值 结果 , 虚线 是 对 应 EPRT (Equation of Phonon Radiative 
Transfer) 的 DOM (Discrete Ordinates Method) 数值 
解 (231, (b) 薄膜 法 向 有 效 热 导 率 随 Kn 的 变化 规律 ,离散 的 空 
心 方 格 点 其 当前 的 MC 数值 结果 , 实 线 是 理论 表达 式 P3] 


Fig. FON Its of cross-plane phonon transport through Si 
K m at room temperature: (a)non-dimensional 


rature distribution at different Kn, symbols are the 


Xo MC results, whereas dashed lines are DOM (Discrete 


ordinates method) solution of EPRT (Equation of phonon 
radiative transfer); (b) effective cross-plane thermal 
conductivity of thin film, squares are the present MC results, 
and solid line is theoretical result[23] 


上 述 结果 表明 ， 本 文采 用 的 动力 型 MC 数值 程 
序 可 以 有 效 捕捉 法 向 和 面向 声 子 输 运 的 非 平衡 效应 ， 
因而 可 用 于 本 文 的 双 层 薄膜 界面 声 子 输 运 研究 ， 


a UR MRO RIRA 
型 并 作 Debye 近似 . 考虑 nails 0.1, 0.4, 1 和 


3 的 四 种 情形 , 统计 宏观 量 的 时 间 间隔 均 为 5x 10-2 
s, 网 格 尺 寸 依次 为 4179x10-9 m、4.179x10-9 
m, 2.0895 1079 m 和 6.9650x10-10 m, 偏差 能 量 粒 
子 数目 分 别 为 3 亿 、3 亿 、4 亿 和 4 亿 . SPER 
截面 均匀 设置 100、 25, 20 和 20 个 网 格 进行 宏观 量 。 CN 
的 统计 ,MC 模拟 结果 如 图 3 所 示 . 结果 表明 , RE \\ 各 0 


热流 分 布 偏离 Fourier 定律 描述 的 均匀 分 布 ， ge Sjal ee 
近 的 热流 受到 抑制 ，Kn 越 大 ， 偏 离 Fourier REE ó y 

程度 越 大 ， 当 前 数值 结果 与 经 典 Fuchs-Sondheimer 02 
模型 的 解析 解 P 吻合 很 好 . 
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图 3 室温 下 硅 薄 膜 中 的 稳 态 面向 声 子 输 运 的 无 量 纲 热流 分 布 示 
意图 : 符号 代表 当前 的 MC 模拟 结果 , 实 线 为 
Fuchs-Sondheimer 模型 的 解析 解 [25] 


Fig. 3 Non-dimensional heat flux distribution of in-plane 


9 MC Kn=0.035 

5 MC Kn=0.35 

* MC Kn=3.5 
=~ Majumdar, 1993-EPRT 


phonon transport through Si thin film at room temperature: 
symbols are present MC results, whereas solid line is the 
theoretical results of Fuchs- Sondheimer model 


3 结果 和 讨论 


3.1 双 层 薄膜 界面 法 向 声 子 输 运 
对 于 双 层 薄膜 材料 界面 法 向 的 声 子 输 运 ,Si/Al 


4 期 FL eS: 双 层 薄膜 界面 声 子 输 运 的 Monte Carlo 模拟 


薄膜 的 总 厚度 为 40 nm， 在 法 向 均匀 设置 40 个 
网 格 进行 宏观 量 的 统计 ， 统 计 宏 观 量 的 时 间 间 隔 为 
5x10-13 s, 网 格 尺 寸 为 1x10-9 m, 式 (7) 中 区 4=300 
K, 模拟 选取 的 总 计算 时 间 为 1.2x107° s, 使 用 的 偏 
差 能 量 粒子 数 为 9000 万 . MC 模拟 的 温度 分 布 结果 
如 图 4 所 示 . 结果 表明 , 在 每 种 材料 内 , 温度 分 布 近 
似 为 线性 , HE Si/Al 界面 处 存在 温度 跳跃 AT, 基 
于 界面 热 导 的 定义 时; 

Gx = 45 (9) 
可 计算 获得 界面 热 导 值 为 Gx=4.2268x10° 
W/(m?K), 和 室温 下 Si/Al 界面 热 导 的 实验 值 
3.45459x108 W/(m? K) 接近 。 这 表明 考虑 真实 
色散 关系 和 频谱 界面 穿 透 系数 的 DMM 模型 能 够 较 
准确 地 描述 界面 声 子 输 运 ,也 验证 了 新 的 MC 数值 
模型 的 可 靠 性 . 
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图 4 室温 下 Si/Al 薄膜 界面 法 向 声 子 输 运 温度 分 布 图 : 圆圈 代 
表 考 虑 真实 色散 关系 和 频谱 穿 透 系数 的 MC 计算 结果 , Si/Al 
薄膜 的 总 厚度 为 40 nm 
Fig. 4 Temperature distribution of cross-plane interfacial 

phonon transport through Si/Al thin film at rod 
temperature: circles are the MC results with exact phonon 
dispersion and frequency-dependent transmission coefficient. 
The thickness of Si/Al thin film is 40 nm 


3.2 双 层 薄膜 界面 面向 声 子 输 运 

对 于 双 层 薄膜 材料 界面 面向 的 声 子 热 输 
运 ，Si/Al 注 膜 总 厚度 为 40 nm， 在 纵向 均匀 设置 
40 个 网 格 进行 宏观 量 的 统计 , 统计 宏观 量 的 时 间 间 
hay 3x107! s, RRA 1x107? m, 式 (7) 中 
Tea=300 K, 模拟 选 取 的 时 间 上 限 为 3x107 s, 使 
用 的 偏差 能 量 粒子 数 是 8000 H. MC 模拟 的 截面 热 
流 分 布 结果 如 图 5 所 示 . 结果 表明 , 每 种 材料 中 的 截 
面 热 流 分 布 类 似 于 单 层 薄膜 中 面向 声 子 输 运 的 截面 
热流 分 布 ， 界面 处 热流 比 薄 膜 内 部 的 热流 更 低 ， 这 
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表明 界面 对 声 子 输 运 的 抑制 作用 与 边界 效果 类 似 。 
另外 , 可 以 发 现 , 截面 热流 分 布 在 界面 处 是 连续 的 . 
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图 5 室温 下 ve 薄膜 中 界面 面向 声 子 输 运 的 截面 热流 分 布 示 
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意图 : 虑 真实 色散 关系 和 频谱 穿 透 系数 的 MC 计算 
R 结果 , Si/Al 薄膜 的 总 厚度 为 40 nm 


eat flux distribution of in-plane interfacial phonon 


he through Si/Al thin film: circles are the MC results 


ith the exact phonon dispersion and frequency-dependent 
transmission coefficients. The thickness of Si/Al thin film is 
40 nm 
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本 文采 用 动力 型 声 子 Monte Carlo 方法 并 引入 
考虑 真实 色散 关系 的 DMM 界面 模型 ， 建 立 了 一 个 


频谱 穿 透 系数 的 使 容声 子 输 运 介 观 数值 模型 对 室 
温 下 Si/Al 褒 层 薄膜 法 向 和 面向 的 界面 声 子 输 运 进 


完 久 弄 面 法 向 声 子 输 运 中 计算 获得 的 界面 热 
结果 吻合 良好 ,新 的 数值 模型 能 准确 捕捉 


行 了 研 


偷 运 的 界面 温度 跳跃 和 面向 输 运 的 非 均 匀 热 流 


At. 本 文 工 作 可 以 很 容易 拓展 至 多 层 薄 膜 的 界面 


声 子 热 输 运 研究 中 。 研究 结果 将 促进 对 界面 声 子 输 
运 机 理 的 认识 . 
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